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In this paper we investigate the electron-phonon contribution to the resistivity of suspended single 
layer graphene. In-plane as well as flexural phonons are addressed in different temperature regimes. 
We focus on the intrinsic electron-phonon coupling due to the interaction of electrons with elastic 
deformations in the graphene membrane. The competition between screened deformation potential 
vs fictitious gauge field coupling is discussed, together with the role of tension in the suspended flake. 
In the absence of tension, flexural phonons dominate the phonon contribution to the resistivity at any 
temperature T with a T 5 ^ 2 and T 2 dependence at low and high temperatures, respectively. Sample- 
specific tension suppresses the contribution due to flexural phonons, yielding a linear temperature 
dependence due to in-plane modes. We compare our results with recent experiments. 
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I. INTRODUCTION 

In recent years the discovery of graphene a mono- 
layer of carbon atoms arranged in a honeycomb lattice, 
stimulated an unprecedented interest from physicists be- 
longing to different communities. Electrons in graphene 
behave as massless Dirac fermions, as confirmed by their 
peculiar Quantum Hall Effect In parallel, graphene 
represents the only existing two-dimensional (2D) con- 
ducting membrane embedded in 3D space. This feature 
became most prominent with the e xpe rimental realisa- 
tion of suspended graphene devicesP^ The latter allow 
for the investigation of the intrinsic properties of the ma- 
terial, unperturbed by the presence of a substrate. Thus, 
suspended graphene offers the unique possibility of ex- 
ploring a system involving at once features of quantum- 
electrodynamics as well as of hard and soft condensed 
matter physicsP 

For a long time the very exi stence of 2D membranes 
was thought to be impossible^ due to their tendency to- 
wards spontaneous crumpling. Indeed, at harmonic level, 
the elastic theory of flat 2D membranes yields divergent 
fluctuations of the angles of the normal vectors to the 
membrane, due to thermal fluctuations of flexural (out- 
of-plane) phonons. It was later understood that the non- 
linear coupling between stretching and bending energies 
hardens the bending stiffness at long wavelengths and 
stabilises the flat phasePEfll 

It was soon realised that mechanical deformations 
of graphene sheets affect the electronic properties 
by inducing coexistent scalar and (fictitious or syn- 
thetic) gauge fields in the effective low-energy Dirac 
HamiltonianP^K^ Recently this issue became of special 
relevance and stimulated the emergence of the so-called 
strain-engineering community, aimed at controlling the 
electronic proper ties o f graphene by suitably engineering 
the deformations.^^ 

Suspended membranes show flexural deformations 



which are typically extremely soft compared to in-plane 
ones. It may be tempting to conclude that the for- 
mer should dominate the low-energy electromechanical 
properties of suspended graphene due to their large den- 
sity of states. This effect is balanced by the strength 
of the electron-phonon coupling in graphene. While in- 
plane phonons show a conventional linear coupling, the 
intrinsic ele ctron- phonon coupling of the flexural modes 
is quadratic! 10 1 12 1 This feature is protected by the re- 
flection symmetry with respect to the plane, as the ef- 
fect of out-of-plane modes cannot depend on the sign 
of the deformations. Thus, we find an interesting com- 
petition between hard-to-excite but strongly-coupled in- 
plane phonons and soft but weakly-coupled flexural ones. 

From an experimental point of view, this competition 
can be addressed, e.g., by transport measurements an- 
alyzing the temperature dependent resistivity Ap. In- 
deed, at high enough electron concentrations and for 
not too low temperatures (where quantum effects be- 
come prominent 1 -), Ap is ascribed to electron-phonon 
scattering. Recent measurements^^ show a Ap scal- 
ing linearly with temperature T. This behaviour is com- 
patible with the expectations for longitudinal in-plane 
phononsP12U a t T > (with T$, the related Bloch- 
Gruneisen temperature) as if the flexural phonon con- 
tribution were irrelevant. In addition, non-suspended 
samples^ show a density independent Ap while sus- 
pended graphene^S exhibit a density dependence which 
has not been addressed so far. 

In a previous paper- we investigated the low- 
temperature dependence of the resistivity of suspended 
graphene in the absence of tension and found that flex- 
ural modes should actually dominate over the in-plane 
ones yielding a T 5 / 2 log T dependence. Motivated by 
the experimental puzzles above, here we further anal- 
yse the different contributions of in-plane vs flexural 
phonons to the temperature dependent resistivity by ex- 
ploring the high-temperature regime and addressing fur- 
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ther issues related to non-universal factors, like contact- 
induced tension in the membrane. In particular we in- 
vestigate whether flexural modes keep dominating the In- 
dependence of the resistivity also at high temperatures 
and how tension affects their contribution with respect 
to in-plane phonons. We discuss the dependence of the 
resistivity on temperature as well as on electron density 
in light of the recent experimental findings. 

The structure of the paper is as follows: In Sec. [TT] 
we introduce the basic description of the electronic and 
phononic properties of graphene to be employed in the 
rest of the manuscript. In Sec. |III| we discuss the in-plane 
vs flexural phonon contributions to the resistivity within 
a Boltzmann approach. We deduce the temperature and 
density dependence of the resistivity Ap in the different 
regimes where acoustic phonons are relevant. We also 
comment on the effect of screening in the density de- 
pendence of the electron-phonon resistivity. Finally, we 
conclude in Sec. IIVI 



II. 



ELECTRONS AND PHONONS IN 
GRAPHENE 

A. Electronic properties 



In this paper we consider the scattering between elec- 
trons in graphene and long-wavelength phonon modes. 
We can thus treat the two Dirac cones independentl y and 
focus on the effective low energy Dirac Hamiltoniar i 22 * 23 * 



H = fiv a • k 



(1) 



where v ~ 10 6 m-s _1 denotes the Fermi velocity and 
the 2D wavevector k is measured from the relevant 
Dirac point. The Hamiltonian in Eq. ([I]) acts on two- 
component spinors (u A k , u B k ) of Bloch amplitudes in 
the space spanned by the two inequivalent honeycomb 
sublattices (A — B). The components of the vector <x 
are the Pauli matrices in the sublattice space. The 
real spin degree of freedom does not play a role here 
and will be ignored. The electronic spinor eigenstates 
|k) with chirality s = er • k/|k| = ±1 have energy 
e k = s?ii> |k|. Their real space representation is (r|k) = 
l/v2(l, s exp[z</> k ]) exp[ik • r], with r the position vector 
in the 2D plane and <j> k the angle of the vector k with 
respect to the x-axis fixed along the armchair direction 
of the graphene lattice. 



B. Phononic properties 

The mechanical distortions of graphene are described 
by the vector field u(r) and by the scalar field h(r) asso- 
ciated with in-plane and flexural (out-of-plane) deforma- 
tions, respectively. The physics of the mechanical distor- 
tions is captured by the elastic Lagrangian density 
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with contributions coming from stretching and bending 
energies. Here po is the mass density of graphene and 
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the strain tensor) 2 ^ with i,j 6 {x, y}. The Lame coef- 
ficients A and p characterise the in-plane rigidity of the 
lattice, while k is the bending stiffness and 7 is a sample- 
specific coefficient associated with tension induced by 
the edges of the sample.^ Notice that the 7(V/i) 2 term 
breaks rotational symmetry which would be obeyed in 
the absence of tension. 

In the harmonic approximation the Lagrangian above 
yields two in-plane phonon modes (longitudinal I and 
transverse i) and one flexural branch (h) with dispersions 
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and group velocities = [(2/x + A) / po] — 2 • 10 4 m/s 
and = [p/pg] 1 ^ 2 ~ 1.3- 10 4 m/s. At long wavelengths 
(with respect to a — 1.42 A, the graphene lattice spac- 
ing), in-plane phonons have a linear dispersion as Gold- 
stone modes associated with the breaking of translational 
invariance in the plane. In contrast, flexural modes have 
a quadratic dispersion in the absence of tension. Ten- 
sion introduces a new wavevector scale = (7/K) 1 / 2 
discriminating a tension-induced linear dispersion at low 
momenta from the quadratic dependence at q > 
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with a = (7/P0) 1 / 2 and f3 = (k/pq) 1 / 2 . Correspondingly, 
the density of states (DOS) for in-plane modes is linear in 
energy e while that for flexural phonons is linear in energy 
up to hujq^ and independent of energy for e > Tiujq^ . 
As we will show, the DOS reduction induced by tension 
is the basic mechanism suppressing the contribution of 
flexural phonons relative to the tension-independent in- 
plane ones. 

In the following it will be useful to Fourier transform 
the phononic displacements x^\r) (with v = l,t,h) as 
x^ u \v) — ^ q Xq exp[iq • r]. The normal modes can 
then be expressed as Xq^ = £q (aq + a^) with = 

[h/2Mujq" > ] 1 / 2 the oscillator length, aq the annihilation 
operator for the mode v at wavenumber q and M the 
total oscillator mass per unit area. 
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C. Electron-phonon coupling 

In this paper we focus on intrinsic coupling mecha- 
nisms between electrons and phonons due to the effect 
of deformations on the electronic Hamiltonian. Other 
sample-specific mechanisms exist, e.g. due to the ca- 
pacitive coupling between a back-gate and electrons in 
suspended graphene or due to buckling of the membrane 
(which breaks the reflection symmetry of the membrane) , 
but these are beyond the scope of the present work. 

The intrinsic coupling between electrons and deforma- 
tions is related to the variation of areas and lengths in- 
duced in the membrane by specific distortions. As the 
variations of length or area are described by the com- 
ponents of the strain tensorp^l by examining its form in 
([3]) one readily concludes that in-plane phonons have a 
linear coupling to electrons while flexural phonons have 
a quadratic one, as long as the reflection symmetry with 
respect to the plane is not broken. We point out that the 
breaking of reflection symmetry (e.g. via capacitive cou- 
pling or via buckling) would yield a non-universal linear 
coupling for flexural modes. In the presence of tension 
this would effectively result in a sample-specific contri- 
bution to the resistivity with identical parametric depen- 
dences as for in-plane phonons. 

In more detail, the representation of the electron- 
phonon coupling in the electronic Dirac description of 
a single valley is given by^ 1 1 
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with /[ity] = 2u xy + i(u xx — Uyy). The diagonal part 
of the coupling constitutes a scalar deformation poten- 
tial originating from local area variations. The corre- 
sponding bare coupling constant has been estimated to 
be^l g 1 ~ 20 — 30 eV. In addition, distortions which in- 
duce no variation of local areas (e.g. pure shear modes) 
would still couple to electrons via the induced bond- 
length modulations. These affect the hopping amplitudes 
between neighbouring carbon atoms and induce the off- 
diagonal terms of Eq. ^ corresponding to a fictitious 
(or synthetic) gauge field in the Dirac equation. The 
corresponding coupling constant has been estimated to 
btP g 2 ~ 1.5 eV which is about an order of magnitude 
weaker than the deformation potential coupling. In the 
following, we will assume that these estimates are at least 
roughly correct. However, we emphasize that it would be 
straightforward to adapt our results to situations where 
this inequality no longer holds. 

It should be noticed however that the deformation po- 
tential, contrary to thegauge field coupling, is affected by 
electronic screening } 26 * 27 ! which reduces the coupling con- 
stant g 1 . Indeed, if Q is the wavevector transferred in 
the electron-phonon coupling, Thomas-Fermi screening 
of the deformation potential yields the coupling constant 
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with Q — |Q|, Vq = 2ne 2 /Q the 2D Coulomb interac- 
tion and IIq the fermionic polarization. The Thomas- 
Fermi screening wavevector Q TF = 27re 2 ^ F = rja^k-p is 
expressed via the fine structure constant of graphene 
a { = e 2 /Tiv ~ 2 and the electronic DOS at the Fermi 
level v F = rjk-p I 2ttTiv , with 77 = 4 due to the spin and 
valley degeneracy. For small wavevectors Q <C Qtf 
the deformation potential coupling is thus strongly sup- 
pressed. We can identify a crossover wavevector Q GD = 
Qtv9iI (9i ~ 92) ^ 2fc F below which the gauge field cou- 
pling dominates over the deformation potential. 

In summary, the longitudinal in-plane and flexural 
modes, as they induce local variations of area, cou- 
ple to electrons via both the deformation potential and 
the gauge field mechanisms, while transverse in-plane 
phonons involve pure shear and couple only via the gauge 
field. At wavevector Q, these phonon modes are charac- 
terized by the coupling matrices 
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in the Dirac description. For flexural modes Q = qi + 
q2, with q! and q2 the wavevectors of the two flexural 
phonons. In the matrices above $, </>i and <f>2 are the 
angles of the vectors Q, qi and q2 with respect to the 
x-axis, respectively, and we defined 4> = 4>\ — 4>2 and 

i> = 4>\ + fo- 



wl. CONTRIBUTION TO THE RESISTIVITY 
FROM ELECTRON-PHONON SCATTERING 

We now address the electron-phonon contribution to 
the resistivity of suspended graphene. It is believed that 
this yields the dominant temperature dependence of the 
resistivity for the reg ime o f high temperatures involved 
in recent experiments^ 1121 (50 K < T < 300 K). Fur- 
ther mechanisms contributing to the temperature depen- 
dence of the resistivity discussed in the literature focus 
on the ballistic regimep^ on the T-dependent screen- 
ing of Coulomb impurities p2l or on the combined ef- 
fects of electron-electron interactions and atomic-scale 



impurities 
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We focus on the high electron density regime where the 
Fermi wavevector fc F is larger than the inverse mean free 
path due to disorder and electron-phonon scattering. In 
this regime a quasiclassical Boltzmann approach to trans- 
port can be employed. The linearized Boltzmann equa- 
tion for the electronic distribution function f k in presence 
of a constant electric field E has the standard fornPS 



(9) 



with — e the electron charge and v = <9e k /<9(?ik) the elec- 
tron velocity. Here we assumed the stationary condition 
df k /dt = as well as spatial uniformity, df k /dr = 0, 
and we denoted by / k = (exp[(e k — n)/k B T] + the 
equilibrium Fermi-Dirac distribution function, with \i the 
chemical potential. The collision integral on the right 
hand side of Eq. ^ depends linearly on the deviation 
from the equilibrium distribution, expressed through the 
function ip k as 



f - f - d/k 
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In the relaxation time approximation we can express the 
collision integral as 



k - /k = ^ M 
r k 9e k 



(ii) 



in terms of the transport scattering time r k . By direct 
comparison of Eqs. ^ and (111 we obtain the solution 
of the Boltzmann equation as 



-eE ■ vr k oc cos£ 
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with £ the angle between E and v. Notice that in 
graphene the modulus of the velocity is independent of 
wavevector and v||k. The solution to the linearized 
Boltzmann equation yields the current density 
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In the regime k B T <C e F (with e F = hvk F the Fermi 
energy), df k /de k is sharply peaked around the Fermi 
level and the current density can thus be obtained by 
performing the angular average of over k with the scat- 
tering time r k evaluated at the Fermi level. As a result, 
j = e 2 v 2 v ¥ T k E/2, yielding the longitudinal resistivity 
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with n = rjkp /Air the electronic density. 

We now focus on the collision integral C(ip k ) for the 
scattering of electrons on in-plane and flexural phonons 
in graphene. By comparison with Eq. (11) we will thus 



deduce the respective transport scattering time, leading 
to the resistivity via Eq. ( 14 ) . 

A. In-plane phonons 

Let us first consider scattering between electrons and 
in-plane phonons of type v (y = I, t for longitudinal and 
transverse modes, respectively). The collision integral 
CW(/ k ) describing the detailed balance of the occupa- 
tion of an electronic state with wavevector k is given by 
the Fermi golden rule expression 



2tt 
T 

2tt 



E W k',Q;k f/k'(l - h)n% ] - / k (l - 40(1 + n<g)] S(e k - e k , - 



Q 



E^U /k<(l-/k)(l + ^V/k(Wk<K} S(e k -e k , + fi^), 



(15) 



where is the distribution function for phonons and 



the factors W k ^. k q are given by 
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The first term describes the process of absorption of a 
phonon with wavevector Q and energy huj^ by an elec- 
tron with wavevector k' and the reverse process involving 
the emission of the phonon by an electron with wavevec- 
tor k. In parallel, the second term describes the emission 
of the phonon by the electron in state |k') and the ab- 
sorption of the phonon by the electron in |k). 



We now proceed to lin eariz e the collision integral 
abovepil making use of Eq. (fl"o|. In doing so, we expand 



e k ,±huj Q 



valid if 



TiWq^ e k . Indeed, although, strictly speaking, phonon 
scattering is inelastic, acoustic phonons in graphene have 
a low group velocity with respect to the electrons, which 
justifies a quasielastic approximation for the scattering 
rate. We also assume the phonons to always remai n in 
equilibrium!^ Summing up the contributions in Eq. ( 15 ) 
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we obtain the linearized collision integral 
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where JIq = (exp[?iWq'/fc B T'] — 1) 1 is the equilibrium 
Bose- Einstein distribution. As the standard solution ( 12 ) 
yields tp k oc cos C, we can write 
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with 9 the angle between the electronic wavevectors k and 
k + Q. By direct comparison with Eq. (Ill, we deduce 
the transport scattering rate 
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The derivative of the Bose distribution implies that the 
relevant phonons to be considered have energies up to 
hujQ^ r~j k B T and their wavenumbers restricted to Q < 
qip\ with TujJ (2) = fcgT". This can be understood as 

follows. In absorption processes, it is only these phonons 
that have a large enough equilibrium occupation, while 
in emission processes, electrons can transfer at most an 
energy of order k B T to the phonons due to the Fermi 
distribution. 



In addition, the quasielastic approximation in Eq. ( 19 ) 



demands energy conservation involving only the incom- 
ing and outgoing electron energies at the Fermi level. 
When combined with the single-valley approximation, 
the quasielasticity implies conservation of chirality. In 
the following, we thus choose s = 1 without loss of gen- 
erality. The on-shell condition and the relation between 
momenta and their angles (sketched in Fig. [TJ yield 
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where <&o denotes the angles between k and Q fulfill- 
ing cos<i>o = — sin(0/2) = —Q/2k F . This condition also 
forces Q € [0, 2k F }. As a result we deduce 



l-cos(9 = 2sin 2 (6>/2) = 2 
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FIG. 1. The on-shell condition forces the incoming and out- 
going electron momenta k and k + Q on the Fermi circle, with 
Q the phonon- induced scattering wavevector. 



graphene, for two reasons: The first reason is conven- 
tional and is associated with the Bloch-Gruneisen tem- 
perature T^q = fu^2k I ■ When T > , large angle 

scattering is possible while, at T -C T B( l, onr y small-angle 
scattering is relevant. The second reason is specific to 
graphene and is associated with the screening dependence 
of the relevant electron-phonon coupling mechanism. If 

T ^> TgQ, screening hardly affects the dominant defor- 
mation potential contribution to the resistivity, while, at 
small enough temperatures, it suppresses the deforma- 
tion potential in favour of the gauge field coupling. 

Indeed, in the high temperature regime T > T$ the 
sum over phonon momenta is cut off at 2k F and scattering 
at all angles is possible. Consequently, screening does not 
significantly affect the deformation potential coupling. 
The latter yields the dominant contribution to the resis- 
tivity for longitudinal phonons with g± (Q) ~ g 1 while 
transverse modes depend only on the unscreened gauge 
field coupling. 

In contrast, in the regime T -C T B( l, the sum is cut off 

at q^ (with q^) <C 2k F ). This constraint implies domi- 
nant small angle scattering (9 <C 1), and the term 1— cos 9 



therefore yields the factor 2{q^ ) /2k F ) 2 ~ {T/T^f < 1. 
In addition, in this regime the deformation potential cou- 
pling for longitudinal phonons is strongly screened. We 
thus identify an intermediate regime for Tq B < T < T b q 
(with = TuJq /k B ) where small angle scattering is 
accompanied by a dominant deformation potential. Fi- 
nally, at low-temperatures T < Tqq the scattering rate 
is dominated by the unscreened gauge field coupling. 
As screening does not affect the coupling for transverse 
phonons, the latter do not show an intermediate temper- 
ature regime. 
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1. Longitudinal phonons 



We thus conclude that the phonon-induced resistivity 
exhibits a number of temperature regimes in suspended 



We first turn to the transport scattering rate due 
to longitudinal in-plane phonons which couple to the 
carriers through both, a gauge field and the deforma- 
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tion potential. The deformation potential contribu- 
tion to th e scatte ring rate has been analysed before for 
monolayer^ 21 * 32 * 33 ! and bilayer grapheneP^l Note however, 
that while the deformation potential is believed to have 
a larger bare coupling constant, it is suppressed at long 
wavelengths due to screening for nonzero doping. Thus, 
screening leads to temperature regimes in addition to 
those associated with the Bloch-Griineisen temperature. 
With the coupling matrix in Eq. (Tsl) we deduce 



+ Q|™g|k> 
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including both deformation potential and gauge field cou- 
plings. The dominant contributions to the scattering rate 
(191 in the various temperature regimes therefore are 
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with T = g\k\ /2hp vvV\ T' = Tg 2 /g 2 and f{x) a func- 
tion fulfilling f(x < 1) ~ 1 - 32x/37r and f (x > 1) ~ 
l/2x 2 . 

The T 4 dependence for T <C Tq^ is due to the gauge 
field coupling. It is analogous to the T 5 dependence 

(Bloch law) for T < in three dimensions pS the dif- 
ference in the power law stemming from the reduced di- 
mensionality of momentum space. 
In the intermediate temperature regime, 

T GD < T < 

TgQ, the screened deformation potential dominates over 
the gauge field coupling. The T 6 dependence originates 
from the combined effects of screening and small angle 
scattering. 

Finally, the linear temperature dependence for T T^q 
stems from the high-temperature expansion of the Bose 
distribution. Since the typical momentum transfers in 
this regime are of order 2k ¥ , it is associated with the es- 
sentially unscreened deformation potential. This regime 
turns out to be the mo st rel evant for the interpretation 
of recent measurements JSEU 



2. Transverse phonons 

The analysis above can be extended to transverse in- 
plane phonons. These have a slightly smaller group ve- 
locity than longitudinal modes and couple only via the 
gauge field mechanism, so that screening is not relevant. 
As a result, the intermediate regime in Eq. (22 1 is absent. 



The coupling term in Eq. ^ yields 
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leading to the scattering rate 
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for T » Ti% (24) 
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with T G = g^kp/hpovv^. The transverse in-plane 
phonons thus yield a contribution to the resistivity com- 
parable to the longitudinal ones for T < T^, while for 
higher temperatures they can be neglected to a good ap- 
proximation. 



3. Temperature- dependent resistivity in recent experiments 

Recent experiments in graphene highlighted the con- 
tribution to the resistivity due to electron-phonon scat- 
tering A/9 at T > TgQ plUH In this regime, using our 
estimates above, the resistivity due to scattering by in- 
plane phonons, Ap in = e F /ne 2 v 2 T k , is dominated by the 
deformation potential coupling of longitudinal modes and 
is given by 



^/(Qtf/^f) 
Ahp e 2 v 2 vW 2 



k B T 



2-10- 8 glf (25) 



in terms of the rescaled quantities g 1 = 
g J(Q TF /2k F )/eV and f = T/K. The contribu- 
tion to the resistivity due to in-plane phonons is thus 
linear in temperature and independent of the electron 



concentration for T T, 



(0 



BG" 



Experiments in non-suspended samples 18 show a 
temperature-dependent component of the resistivity 
compatible with Eq. (251, which is quantitatively consis- 



tent with the estimates for the bare deformation poten- 
tial coupling constant.^ Suspended graphene devices^ 
likewise show a linear T dependence, compatible with 
in-plane phonons, but accompanied by an unexpected 
electron-density dependence not captured by our analy- 
sis above. In particular, for increasing electron concen- 
tration the resistivity has been observed to decrease and 
finally saturate to the in-plane phonon contribution ( 25 ) 
at large densities. 

While it may be tempting to think that this behaviour 
originates from the screening of the deformation potential 
at increasing electron density, our analysis shows that 
this is not the case. Indeed, for T > T^, Q is cut off at 
2fc F . This, together with the fact that <5 TF oc fc F , implies 
that the screening factor in Q does not introduce any 
additional density-dependence into the scattering rate. 
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Moreover, if this was a relevant effect, it should appear 
in non-suspended samples as well. 



The in-plane phonon contribution ( 25 ) seems suffi- 



cient to describe the behaviour of non-suspended sam- 
ples, where flexural (out-of-plane) deformations are sup- 
pressed by the direct contact with a substrate. In 
graphene samples on very rough substrates small out- 
of-plane fluctuations could still survive due to the non- 
perfect adhesion of the membrane to the surface of the 
substrate. 

In contrast, flexural phonons can become relevant in 
suspended devices. We now turn to analyse their effect 
along the same line as for the in-plane modes above. 



B. Flexural phonons 

The analysis of the resistivity due to flexural phonon 
modes is more involved due to the presence of two 
phonons with wavevectors q 1 and q 2 at each interac- 



tion vertex. These yield four possible processes involving 
double absorption, double emission or mixed absorption- 
emission terms in the scattering rate. In addition, four 
relevant wavevector scales come into play, namely q^\ 
Q GD , and hp, related to temperature (turfy^ = k B T), 

screening, tension, and electron concentration, respec- 
tively. A further ultraviolet scale, corresponding to a 
Debye energy of order hcu^ a ~ 100 meV, is assumed to 
provide the largest energy scale. For temperatures higher 
than hojy^/k-Q the straightforward expansion of the two 
Bose distributions yields a scattering rate depending on 
temperature as T 2 . However, at those energies the elastic 
treatment of the membrane as well as the elastic approx- 
imation in the scattering rate may be questionable. We 
thus focus on lower temperatures which are relevant to 
experiments. 

The collision integral including all emission and ab- 
sorption processes for scattering of electrons and flexural 
phonons has the form 



cih) (h) = y E w Sl,^ M 1 - /kK'X^ - /k(i - /*)(! + + <>) I - - H h ? - M?) + 



qi.q 2 



2tt 
2n 



E <tq 2 ;k< [/k'(l - /k)(l + + <) - 4(1 - fv)ngng] S(e k - e k , + hug + Hug) + 



E < } 

qi>q 2 



k^k'^ 



4<(1 - 4)<(1 + < } ) - 4(1 - 40(1 + ng)ng] S (e 



E W k,q i; k',q 2 4'(1 - 4)(1 + - 4(1 - 40< ) (1 + <) S(e k - e k , + hug - hug) (26) 



qi>q 2 



with 



W, 



kik'.q^q.-, "k.k'+qj^+qj. 



(k|<! q Jk') 



(27) 



A lengthy but straightforward calculation analogous to 
the case of in-plane modes above yields the linearized 



collision integral of the form (II), with the transport 
scattering rate 



1 



2tt 



dug <9w q . 



q^q 2 



CO 



h 



(1 - COS0) X 



Q|<i q Jk)|*5(e k+Q -e k )x 



(28) 



ui, 



(h) 

0(h) 



u, 



(h) 



0(h) 



.M 
^q t 

0(h) 



Xh) 
0(h) 

Mi J 



with 9 the angle between the electronic wavevectors k and 
k + Q and Q = q ± + q 2 the total wavevector transferred 
to the phonons in the scattering process. 

The general considerations in Sec. |III A| hold for flexu- 
ral phonons as well. In particular, the on-shell condition 



( 20 1 as well as Eq. ( 21 1 remain valid. The matrix element 
of the electron-phonon coupling is given from Eq. (Tsl) as 



( k l<!qJ k - 



Q) 



g[ BC \Q)co S <p(l + e ie ) 



,iW>+0 k +e) 



_ e -i(ip+<t> k ) 



Temperature enters the scattering rate ( 19 1 via the 



Bose distributions and cuts off the relevant phonon mo- 
menta at q±, </2 ^ <1t ■ Contrary to what happens for 
in-plane modes, this is true even for qg ^> 2fc F (i.e. for 
T > T^}, with T$ = hug p /k B the Bloch-Griineisen 

temperature of flexural phonons). Indeed, while the to- 
tal scattering wavevector Q is bounded by the on-shell 
condition to fulfil Q < 2fe F , each individual qi and c?2 

can be large, provided that |qi — q 2 | < qg . In this 
regime, scattering at all angles is possible and the domi- 
nant electron-phonon coupling mechanism is via the un- 
screened deformation potential. 

In the opposite limit of low temperature where T <C 



r, 



3Q (i.e. for q^ <C 2k F ) the scattering rate is dominated 

by small angle scattering, with 1 — cos (9 ~ {q^ /2k F ) 2 <C 
1. In this regime, screening sup presses the deforma- 



9iQt /Q 



tion potential coupling in Eq. ([29]), yielding g± (Q) 



TF 



•C <?!■ This suppression favours the un- 
screened gauge field coupling g 2 for -C Q GD (i.e 
for T <C Tqq , with T^d 



^<£_/*b)- 



Tension affects the phonon spectra yielding a linear 
dispersion for q < and a quadratic one for q > q*. 
As the density of phonon states is much higher in the 
region of the quadratic spectrum, the dispersion relation 
which dominantly contributes to the resistivity is linear 
if q^ -C q„ and quadratic for q^ 3> ?*. This DOS 
suppression at low energy is not accompanied by a change 
in the flexural-phonon coupling which remains quadratic 
even in presence of tension as the reflection symmetry 
with respect to the plane is preserved. This is the basic 
mechanism suppressing the contribution to the resistivity 
of flexural-phonons with respect to in-plane ones. 

It is possible to deduce the relevant temperature and 
electron concentration dependence of the scattering rate 
by simple power counting, once the relevant momenta 



are identified. Indeed, for q 



< 



'i B T/huJ { q 



(h) 



(h) 



one has 



Ah) 



and the scattering rate is given by 



qi>q 2 



k\T 2 



h 2 ui^ 'wi' 1 ' 



(1 — cos 9) x 



(29) 



k + Q|«;W \k)\'8(e 



k+Q 



where the factor 4 accounts for all emission-absorption 
processes of the two phonons. Introducing the total and 
relative momenta Q and q = qi — q2, the transport scat- 
tering rate for electrons at the Fermi level scales as 



1 



T 



dqq 



r C« 



.2Av 



dQ x 



(30) 



Q 



Q + Q 



TF 



2 A 



2 2 



, ,(h)2 (h)2 



Here A = for T <C Tqq when the gauge field coupling is 
dominant, while A = 1 for T 3> T G d when the screened 
deformation potential dominates. In principle seven dif- 
ferent regimes can be identified, according to the relative 
magnitude of the relevant wavevectors. 

At low temperatures T <C Tqd ^ ^bg > corresponding 

4 (h) 
to qj, 



<C Qqd ^ 2fc F , A = and small angle scat- 



tering dominates the scattering rate. We can rescale all 
wavenumbers by q^ and find two regimes: 



I: For q, « 
yielding q£ ] 



the relevant dispersion is ui 
VT and l/r k ~T b / 2 /k 2 . 



As pointed out before,^ the first regime with low tension 
shows a logarithmic infrared singularity. For a sufficiently 
large membrane, these singularities are regularised by 
a temperature-dependent infrared cutoff related to the 
non-linear coupling of bend ing and stretching energies in 
the elastic LagrangianP 1 ^ Small tension or finite size ef- 
fects yield alternative infrared cutoffs modifying the value 
of the scattering rates by numerical prefactors of order 
one. 

The expansions above show that, in the first case with 
very weak tension (q^ <C q^ <C Qgd)> the scattering 
rate due to flexural phonons has a T 5 / 2 temperature de- 
pendence and is larger than the corresponding one due 
to in-plane modes, which scales as T . In contrast, for 

q^ « 9, « Qgdi tension suppresses the flexural con- 
tribution to the subdominant T 7 term, which allows in- 
plane scattering to dominate the clcctron-phonon scatter- 
ing contribution to the resistivity at low temperatures. 

At intermediate temperatures T G q «T« ?b G , cor- 
responding to Q GD <C q^ <C 2k F , small angle scattering 
still dominates the scattering rate while A = 1. Screening 
thus yields an additional factor (q^ /Qtf) 2 m the scat- 
tering rate. Rescaling all wavenumbers by q^' we find 
two regimes similar to the previous analysis: 

• III: For q„ <C q^ the relevant dispersion is ui ~ q 2 , 
yielding q% ] ~ y/T and l/r fc ~ T 7 ' 2 /k F . 

• IV: For q)^ ' -C q* the relevant dispersion is ui ~ q, 
yielding q^ ~ T and l/r fe - T 9 /k F . 

Finally, in the high temperature limit T ^> Tg G , corre- 
sponding to q^ ^> 2fc F , large angle scattering is possible 
for qi, qi > fc F , which yield the dominant contribution 
to the scattering rate. In this regime A = 1 and the total 
wavevector Q is cut off at 2k F due to the on-shell con- 
dition while q is limited by q^p ■ In the case of strong 
tension q^ ^> the suppression of the phononic DOS 
at low energy shifts the dominant momenta to large val- 
ues of order q ~ qip. In contrast, small tension <C q^ 
tends to favour low momenta down to q t or fc F , whichever 
is smaller. We identify three regimes: 

• V: For q^ <C 2fc F <C the relevant dispersion is 
us ~ q 2 , yielding a dominant contribution for small 
wavenumbers down to q ~ fc F resulting in l/r fc ~ 



T 2 /k, 



• VI: For 2fc F <C q* <C q^ the relevant dispersion is 
ui ~ q 2 , yielding a dominant contribution for small 
wavenumbers down to q ~ g t resulting in l/r fe ~ 

T 2 fc„. 



II: 



For 



yielding g^ 1 ' 



<C the relevant dispersion is to 
T and l/r t ~T 7 /k 2 . 



q, 



• VII: For 2k F <C q^ <C the relevant dispersion is 



ui 



q, yielding q^ ~ T and l/r fc ~ T 4 k F . 
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FIG. 2. The dependence of the resistivity due to scattering off 
fiexural modes on temperature T and electron density n. The 
grey area identifies the region q„ > where the relevant 
fiexural phonon dispersion is dominated by tension and w, ~ 
aq. 



It has to be noted that the T 2 scaling in regime V stems 
from the relevant q 2 dispersion of fiexural modes with 
low tension, and is not trivially obtained from the high- 
temperature expansion of the two Bose distributions, un- 
like in the case of in-plane phonons. In this case, in fact, 

the qi and (72 momenta are still limited by q^ . 

The corresponding dependence of the resistivity on 
temperature and electron density in the seven regions 
above is summarised in the diagram of Fig. [2] 
By a more explicit evaluation of Eq. (28) we calculate 



the scattering rate due to electron-phonon coupling for 
flexural-modcs in the seven regimes above, yielding 



1 
1 

1 

T k 

F 

1 



(k B T) 



5/2 



(k B T) 7 



C 



G (?k*) 5 a 4 (2fc F ) 2 

(k B T) V2 
{hpf/ 2 ^(2k F Q TF ) 



C 



C 



C 



C 



(ha) 7 a 4 {2k F Q 



(k B Tf 

(k B Tf 2k, 
PHI 

(k B T) 4 2k, 
a 4 (ha) 2 



in region I 



in region II 



in region III 



in region IV 



in region V 



in region VI 



in region VII 



with C ~ gl / (2ir) 3 PqH v and C G ~ Cg\jg\, up to nu- 
merical prefactors of order one. 



C. In-plane vs fiexural phonons 

Our analysis allows us to compare the contribution 
to the resistivity of in-plane and fiexural modes. In 
order to quantify the importance of these two contri- 
butions at a given temperature, it has to be pointed 
out that the Bloch-Griineisen temperatures for in-plane 
and fiexural phonons can be significantly different, in 
particular in the low-tension regime where out-of-plane 
modes have a soft quadratic dispersion. For graphene 
one finds T b L q ~ 50n 1 / 2 K while (in the absence of ten- 
sion) TgQ ~ 0.47T.K, with n = n/10 12 cm~ 2 the rescaled 

electron density. In the absence of tension at T = T b q 
the estimates above yield a ratio between the scattering 
rate due to fiexural and in-plane phonons of 2.5 n -1 / 2 
and for higher temperatures the fiexural modes will be 
even more dominant. 

In practice, at the present electron concentrations for 
suspended graphene samples, in the absence of tension 
the flexural-phonon contribution to the resistivity should 
dominate over the in-plane one at any temperature, show- 
ing a crossover between a T 5 / 2 to a T 7 / 2 dependence 



around Tqd' an< ^ between T 1 ! 2 and T 2 around T^j. 
The main reason why the effect of fiexural modes does 
not appear in experiments is due to the suppression of 
the fiexural phonon contribution to the resistivity by the 
sample-specific tension. Actually, due to the negative 
thermal expansion coefficient of graphenepS tension is it- 
self a temperature dependent quantity. This suppression 
leaves the in-plane contribution as the dominant scatter- 
ing mechanism, yielding a linear-T dependence compati- 
ble with experiments. 

When tension is present, in the experimentally relevant 
regimes q^ <C 2fc F <C q^ and 2k F <C -C the 
contribution to the resistivity from fiexural phonons is 
approximated by 



Ai -9^Lf 2 



n - 



7 



(31) 



with the rescaled tension 7 = 7/2 ■ I0~ 2 Kgs~ 2 . The 
value 7 = 1 corresponds to the condition q^ = 2k F at the 
density h = 1. Comparing Eq. p5|) and (31) we thus get 



Ap 



flex 



1 



T 



Aft, 



20 n + j 



(32) 



We can then estimate the minimal tension needed to sup- 
press the fiexural contribution by imposing the ratio in 
Eq. (32) to be smaller than unity at room temperature. 



Considering that in typical current suspended samples 
h < I, this results in a tension 7 ~ 15, corresponding to 
a strain of about 10~ 3 . These estimates are reliable as 
long as we do not enter the regime of very high tension 
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FIG. 3. The combined contributions to the resistivity due 
to in-plane and flexural-phonons Ap as a function of the 
temperature T for three different electron densities h = 
0.05, 0.15, 0.3 (dashed-dotted, dashed and continuous line, 
respectively). Here we assume a tension 7 = 1 and a de- 
formation potential coupling g 1 = 10. Inset: same plot as in 
the main figure, but for stronger tension 7 = 20. Notice the 
almost perfect linear-T scaling, independent of density. 



2k F <C qj, <C q^, i.e. for 7 < 3T, which is usually easily 



fulfilled for T > T^ G . 

Thus rather weak tension is sufficient to suppress the 
flexural contribution in favour of the in-plane one. This 
is true for temperatures up to T c ~ 20(n + 7) where the 
crossover between the in-plane dominated to the flexural- 
phonon dominated resistivity takes place. The observa- 
tion of the crossover between these two regimes would 
provide information about the otherwise unknown value 
of tension in the sample. This prediction could be tested, 
for example, in graphene samples mounted on break junc- 
tions where tension can be controllably tuned or in flakes 
clamped on a single side with an STM tip as a drain con- 
tact. 

In Fig. [3] we plot the temperature-dependent compo- 
nent of the resistivity Ap = Ap in + Ap flcx according to 
Eqs. ((25| and (31). The combined contributions of in- 



plane and flexural phonons are shown for different den- 
sities in the experimentally relevant range. Even in the 
presence of tension of the order of 7 = 1, the contri- 
bution from flexural-phonons can still significantly affect 
the value of the resistivity. This results in a slight devia- 
tion from the purely linear-T dependence due to in-plane 
modes. In this case, a residual density dependence is ob- 
servable, stemming from the regime q* <C 2fc F <C q^\ 
in qualitative agreement with experiments on suspended 
grapheneP^ In contrast, the inset in Fig. [3] shows the 
density-independent T-linear resistivity at larger tension 
(7 = 20). 

A quantitative understanding of the density depen- 
dence observed in experiments would require the knowl- 
edge of the sample-specific tension, as well as the in- 
clusion of temperature dependent screening of charged 



impurities of Altshuler-Aronov corrections,^ and pos- 
sibly of further non-intrinsic electron-phonon coupling 
mechanisms (e.g. capacitive coupling to a back gate as 
well as buckling). These issues are beyond the scope of 
the present manuscript. 



IV. CONCLUSIONS 

In summary, the interplay between the electronic and 
phononic degrees of freedom in suspended graphene 
membranes offers a rich scenario which can be addressed 
in current transport measurements. 

Here we analysed the contribution to the resistivity of 
suspended graphene due to electron-phonon scattering. 
We discussed the competition between acoustic in-plane 
and flexural distortions in various temperature regimes. 
We focused on the intrinsic electron-phonon coupling in 
graphene due to the interaction of electrons with elastic 
deformations, taking into account both the (screened) 
deformation potential and the fictitious (or synthetic) 
gauge field coupling. Further non-universal coupling 
mechanisms exist, e.g. via the capacitive interaction of 
the membrane with a back gate or the breaking of the 
reflection symmetry (e.g. due to buckling). While we do 
not discuss these in the present work, they yield a sample- 
specific linear coupling for flexural phonons which, in 
presence of tension, would result in a T-dependent con- 
tribution to the resistivity analogous to in-plane modes. 

We find that, for the electron densities achievable in 
suspended graphene, flexural phonons should dominate 
over in-plane ones at any temperature in the absence of 
tension. In particular, the component of the resistiv- 
ity due to flexural phonons should show a T 5 / 2 depen- 
dence for T <C associated with the fictitious gauge 
field coupling. For < T < T$ the T 7 / 2 depen- 
dence stems from the dominant coupling via the screened 
deformation potential while for T T^q screening is 
irrelevant and the resistivity scales as T . A sample- 
specific tension induced by the contacts yields a stiffening 
of the flexural dispersion, corresponding to a suppressed 
phonon density of states, while the reflection symmetry 
with respect to the plane protects their weak quadratic 
coupling. As a result, tension suppresses the flexural- 
phonons contribution to the resistivity. We point out 
that, due to the negative thermal expansion coefficient 
of graphene, tension in real suspended samples is itself 
temperature-dependent. 

We conclude that it is due to the non-universal tension- 
induced suppression of the contribution due to flexural- 
modes that experiments seem to show a resistivity dom- 
inated by in-plane phonons alone. The latter yield a lin- 
ear temperature dependent resistivity for T ^> Tqq, in 
qualitative agreement with experiments. This contribu- 
tion is however independent of electron density even in 
the presence of electronic screening and cannot account 
for the observed density dependence in suspended sam- 
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pies. A density dependence is induced as long as the 
flexural phonons contribute significantly to the resistivity 
which is the case for sufficiently weak tension. It would 
be interesting to probe these issues in samples with a 
controllable degree of tension, like suspended graphene 
in break-junctions or in suspended flakes clamped on a 
single side with an STM-tip as a contact. 



financial support through SFB 658, SPP 1459 of the 
Deutsche Forschungsgemeinschaft as well as DIP. One of 
us (FvO) thanks the Aspen Center for Physics for hospi- 
tality during the final stages of this work. 



ACKNOWLEDGMENTS 



Useful discussions with Eduardo Castro and Maresa 
Rieder are gratefully acknowledged. We acknowledge 



1 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. 
Zhang, S.V. Dubonos, I.V. Grigorieva and A. A. Firsov, 
Science 306, 666 (2004); K.S. Novoselov, A.K. Geim, S.M. 
Morozov, D. Jiang, M.I. Katsnelson, I.V. Grigorieva, S.V. 
Dubonos and A. A. Firsov, Nature 438, 197 (2005). 

2 Y. Zhang, Y.-W. Tan, H.L. Stormer and P. Kim, Nature 
438, 201 (2005). 

3 J.C. Meyer, A.K. Geim, M.I. Katsnelson, K.S. Novoselov, 
T.J. Booth and S. Roth, Nature 446, 60 (2007). 

4 J. Scott Bunch, A.M. van der Zande, S.S. Verbridge, I.W. 
Frank, D.M. Tanenbaum, J.M. Parpia, H.G. Craighead 
and P.L. McEuen, Science 315, 490 (2007). 

5 E-A. Kim and A.H. Castro Neto, Europhys. Lett. 84, 
57007 (2008). 

6 R.E. Peierls, Ann. I. H. Poincare 5, 177 (1935). 

7 L.D. Landau, Phys. Z. Sowjetunion 11, 26 (1937). 

8 D.R. Nelson and L. Peliti, J. Phys. (Paris) 48, 1085 (1987). 

9 M. Paczuski, M. Kardar, and D.R. Nelson, Phys. Rev. Lett. 
60, 2638 (1988). 

10 E. Mariani and F. von Oppen, Phys. Rev. Lett. 100, 
076801 (2008); ibid. 100, 249901(E) (2008). 

11 L.M. Woods and G.D. Mahan, Phys. Rev. B 61, 10651 
(2000). 

12 H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 (2002). 

13 M.A.H. Vozmediano, M.I. Katsnelson and F. Guinea, 
larXiv: 1003.51791 

14 V.M. Pereira and A.H. Castro Neto, Phys. Rev. Lett. 103, 
046801 (2009). 

15 F. Guinea, M.I. Katsnelson and A.K. Geim, Nature 
Physics 6, 30 (2009). 

16 A. A. Kozikov, A.K. Savchenko, B.N. Narozhny and A.V. 
Shytov, i arXiv:1004.0468l 

17 S.V. Morozov, K.S. Novoselov, M.I. Katsnelson, F. 
Schedin, D.C Elias, J. A. Jaszczak and A.K. Geim, Phys. 
Rev. Lett. 100, 016602 (2008). 



18 J.-H. Chen, C. Jang, S. Xiao, M. Ishigami and M.S. Fuhrer, 
Nat. Nanotechnol. 3, 206 (2008). 

19 K.I. Bolotin, K.J. Sikes, J. Hone, H.L. Stormer and P. Kim, 
Phys. Rev. Lett. 101, 096802 (2008). 

20 See, e.g. W.A. Harrison, Solid State Theory, (McGraw-Hill, 
New York, 1970). 

21 E.H. Hwang and S. Das Sarma, Phys. Rev. B 77, 115449 

(2008) . 

22 PR. Wallace, Phys. Rev. 71, 622 (1947). 

23 D.R DiVincenzo and E.J. Mele, Phys. Rev. B 29, 1685 
(1984). 

24 L.D. Landau and E.M. Lifshitz, Theory of Elasticity, 
(Pergamon, New York, 1986). 

25 Typical parameters for graphene are k = 1 eV, fj, ~ 4A = 
9eVA" 2 , p ~ 7.6- 10- 7 kg-m- 2 . 

26 E.H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 
(2007). 

27 F. von Oppen, F. Guinea and E. Mariani, Phys. Rev. B 
80, 075420 (2009). 

28 M. Mueller, M. Brauninger and B. Trauzettel, Phys. Rev. 
Lett. 103, 196801 (2009). 

29 E.H. Hwang and S. Das Sarma, Phys. Rev. B 79, 165404 

(2009) . 

30 V.V. Cheianov and V.I. Fal'ko, Phys. Rev. Lett. 97, 226801 
(2006). 

31 L.D. Landau and E.M. Lifshitz, Physical Kinetics, (Perg- 
amon, New York, 1986). 

32 L. Pietronero, S. Strassler, HR. Zeller and M.J. Rice, 
Phys. Rev. B 22, 904 (1980). 

33 T. Stauber, N.M.R. Peres and F. Guinea, Phys. Rev. B 
76, 205423 (2007). 

34 J.K. Viljas and T.I. Heikkila, Phys. Rev. B 81, 245404 

(2010) . 

35 W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames 
and C.N. Lau, Nat. Nanotechnol. 4, 562 (2009). 



